New set of measures to analyze dynamics of non-equilibrium structures 



^1" 
O 

o 

(N 
5h . 

< 
(N 



43 

o 



CO 

o3 



i 

c 

o 
o 



in 

o 
o 



I 

c 

o 
o 



x 

S3 



Girish NatharQ and Gemunu Gunaratnjj] 

Department of Physics, University of Houston, 
4800 Calhoun Blvd., Houston, TX 77204, USA 

(Dated: February 2, 2008) 

We introduce a set of statistical measures that can be used to quantify non-equilibrium surface 
growth. They are used to deduce new information about spatiotemporal dynamics of model systems 
for spinodal decomposition and surface deposition. Patterns growth in the Cahn-Hilliard Equation 
(used to model spinodal decomposition) are shown to exhibit three distinct stages. Two mod- 
els of surface growth, namely the continuous Kardar-Parisi-Zhang (KPZ) model and the discrete 
Restricted-Solid-On-Solid (RSOS) model are shown to have different saturation exponents. 
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In recent years, there has been a considerable amount 
of effort to analyze non-equilibrium interfacial growth 
and pattern formation in experimental and model sys- 
tems 0, S H H 11 H 0, II H, GjJ. The phenomena 
studied include spinodal decomposition 0, lJEIJEL 
chemical pattern formati on , surface growth 0,13,0)13 
and epitaxial growth [j| llCj . Microscopic modeling of 
these phenomena is highly complex, and most micro- 
scopic details of such structures depend on initial con- 
ditions and stochastic effects. Hence model systems are 
often used to extract statistical properties of these struc- 
tures and to determine the physical processes that are 
most relevant for their growth. Some of the model sys- 
tems that have been introduced to study the spatio- 
temporal dynamics of the aforementioned systems are the 
Cahn-Hilliard Equation (CHE) Q, the Kardar-Parisi- 
Zhang (KPZ) model 0, the Restricted-Solid-On-Solid 
(RSOSVmodel and the Swift-Hohenberg Equation 
(SHE) H, 0. 

In order to confirm if a given model can accurately 
represent a pattern forming process, it is necessary to 
compare as many statistical measures as possible. Such a 
comparison between model systems and physical systems 
is also needed to validate claims of universality. Unfor- 
tunately, there are only a handful of measures that are 
available to be used for such comparisons. The aim of 
this work is to introduce a new family of such measures. 

Commonly used statistical measures to analyze surface 
growth and patterns include surface roughness (i.e., the 
standard deviation of the heights) Wl (t) , where L is the 
lattice size and t the time, the correlation length Q , and 
the domain size [ll|, O For thc KPZ and RSOS inter- 
faces, W L (t) ~ t p when t < L z . At very large times, 
W L {t — ► oo) ~ L a 8]. The growth exponent (3, the dy- 
namic exponent z, and the roughness exponent a, depend 
only on the dimensionality of the growth process and are 
independent of L apart from finite-size scaling corrections 



It has been found numerically that these exponents 
are the same in all dimensions for the KPZ and the RSOS 
models. Based on the results, it has been suggested that 
KPZ and RSOS belong to the same universality class. 
This is a non-trivial statement, since although both mod- 
els consider the competition between random deposition 
and diffusion, RSOS is a discrete counterpart of the KPZ 
with distinct microscopic dynamics. The availability of 
additional statistical measures can be used to test the va- 
lidity of the assertion of universality. In fact, one of our 
conclusions, based on the new measures, is that there 
are statistical features that are not common to KPZ and 
RSOS models. 



Non-equilibrium pattern formation and dynamics has 
also been extensively studied in the context of spinodal 
decomposition using the Cahn-Hilliard Equation 1, 11, 

. The spatio-temporal dynamics can be classified into 
two temporal regimes - an early stage where there are 
many small clusters (cluster size <C L) and a late stage 
where there are large interconnected domains and cluster 
sizes are comparable to L. The domain size in the late 
phase is found to grow in time with an exponent | |l2j| . 



Consider a planar pattern represented by a scalar field 
U(x,t) which can be the height of an interface, or some 
relevant intensity field. One feature not captured by the 
correlation length and roughness is the amount of curva- 
ture of the contour lines of U . The possibility of using 
curvature for such an analysis has been proposed before 
0. E> E E3 » but in practice the results are very sen- 
sitive to noise. The underlying reason is that the evalu- 
ation of k — {U XX U 2 + U yy U 2 - 2U xy U x U y )/ {U 2 X + Uf)i 
is very sensitive to errors in the denominator, when it 
is small. In place of k, one can use another measure, 
namely the determinant of the Hessian normalized by 



the variance of U, A = (U xx U yy ~ U^ y )/Var{U). Un- 
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like k, the calculation of A is fairly insensitive to noise. 
Further, for typical local structures, A is proportional to 
k. For example if U(r) — e~ ar , A = —a 3 e~ 2ar k; for 
U(r) = sin(kr), A = —k 3 sin(2kr)K/2. 
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The new measures we define are, 

^)=(L^J • (1) 

Note that for each 8, (i(8, t) has dimensions of inverse 
length. Furthermore, the measures (i(S, t) are invariant 
under all rigid Euclidean transformations (i.e., transla- 
tions, rotations, and reflections) of a pattern. The use 
of moments, 8, allows us to emphasize regions with dif- 
ferent values of A, thereby providing an array of length- 
scales associated with the structure. The use of multi- 

{ile lengthscales is similar in spirit to using multifractals 
ISL llflj to characterize strange attractors, although the 
spatio-temporal nature of the dynamics for the models 
considered in this Letter makes a further connection dif- 
ficult. In our analysis, we define a(8, t) as the growth 
rate of (i(8,t), i.e., a(8,t) = log((i(8,t))/log(t). 

The organization of the paper is as follows. We first 
discuss our results of our analysis for the CHE, and de- 
scribe distinct stages in the spatio-temporal dynamics 
using (i(8,t). We then proceed to analyze the KPZ and 
RSOS models and show that during the initial stage, 
(i(8,t) does indeed lend additional credence to the sug- 
gestion that in the strong non-linear coupling limit, they 
do belong in the same universality class. However, our 
analysis of the late stage shows that the saturation ex- 
ponents for the two models are different. The computa- 
tional techniques used to obtain the data for these mod- 
els is discussed in detail elsewhere 01 an d wm n °t be 
expanded upon here. 

The CHE models spinodal decomposition using the dy- 
namics of a conservative field ip(x,t) via 

^ = \v 2 (^-^ + ^) + ^li 1 (x,t). (2) 

Here, 77 is delta-correlated noise of zero mean and am- 
plitude 1. The Crank-Nicholson method 20] used to in- 
tegrate Equation (2) allows us to choose timesteps as 
large as 0.4 and enables us to investigate the dynamics 
of the CHE to times large enough to observe saturation 
of fjt(5,t), e.g., T = 800,000 for the lattice size L = 32 
and T = 1,600,000 for the lattice size L = 128. To ob- 
tain (J,(6, t), we used the averages of 32 runs for L = 32, 
16 runs for L = 64, and 8 runs for L — 128. The error- 
bars on all curves are calculated by averaging over these 
realizations. 

The dynamics of domain growth is as follows: Begin- 
ning from a random initial configuration, \ip(x, t)\ and 
the domain size grow in time. For sufficiently lar ge t , 
ip(x,t) reaches its equilibrium values of —1 or +1 |ull2|. 

The behavior of fi(5,t) is shown in Fig. 1(a). The dy- 
namics can be separated into three stages (marked Stage 
1,2, and 3 in the figure). Initially, (i(6, t) appears to grow 
as a power law in time for about one and a half decades, 
as seen in Fig. ^a). There is a formation of small do- 
mains from the random initial field. As seen from the 
topmost curve of Fig. ^b), the growth rates ai(S, t) for 
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FIG. 1: Plots of the measure in time for various moments. 
L — 32. (a) The noise-free case for different 8, (b) The slopes 
versus 8 for Stage 1,2, and 3 for the noise-free case. 



all 8 are nearly identical and independent of t . During 
this stage, one contribution to (1(8, t) comes from points 
near the domain boundaries. This is seen by plotting the 
height along any row or column. Although the formation 
of domains reduces the number of interfacial points , this 
decrease is countered by the fact that since i/j(x,t) has 
not reached its equilibrium values, the internal points of 
a domain provide an increasing contribution to the inte- 
grand. Consequently, (i(8,t) increases during this stage. 

The crossover between Stages 1 and 2 occurs when the 
field ip saturates to its equilibrium values. A histogram 
of interfacial heights clearly shows the concentration near 
the equilibrium values of ±1. Beyond this stage, interior 
points of a domain no longer contribute to (i(8,t). As a 
result, the aforementioned decrease in interfacial bound- 
ary points now leads to a decrease in (i(S,t). The slopes 
(T2 (8, t) in this stage are plotted in the lower curve of 
Fig. Gib). We find that in this stage, the distinct mo- 
ments relax at different rates (i.e, (12(8, t) is 5-dependent). 
Moreover, it is seen from Fig. ^b) that the rate of re- 
laxation of the larger moments is smaller than that for 
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FIG. 2: Bubble disappearance for the L — 128 lattice. The 
actual time is obtained by scaling the x-axis values by a factor 
10 s . Notice how the peak in the measure corresponds to the 
bubble disappearing at approximately time 148, 000. This is 
for 5 = 3.4. 



the lower moments, implying that sharper features (those 
with larger curvature) change relatively slowly. For ex- 
ample, straightening of a domain boundary wall occurs at 
a faster rate than elimination of a small domain. Running 
averages reveal that <7i [8, t) is uniform until the end of 
this stage. Observe one critical difference between Stages 
1 and 2, namely that the difference in relaxation of the 
distinct lengthscales cannot be determined without an 
array of measures. 

The domain size as measured from the two point cor- 
relation function [T^ becomes comparable to L between 
Stages 2 and 3 and 02 (8, t) changes during the transi- 
tion. The spatio-temporal dynamics beyond this involves 
a coarsening of the large domains and elimination of the 
small ones. The disappearance of small domains can be 
identified by peaks in /i(<5, t) for large values of 8. These 
peaks result from the large curvature of a disappearing 
domain. Fig. [21 shows the effect of such an event on 
fi(8,t). It is clear that the elimination of a domain is ac- 
companied by a peak in fi(8,t). The decay rates (7^(8, t) 
for the third stage are not uniform in time, unlike the pre- 
vious stages. An example is shown in the central curve 
of Fig. P^b). Although 03(8, t) varies in time, the general 
form is the same; there is a (5-dependence for small val- 
ues of 8 but a saturation for higher values. At the end of 
Stage 3, only a handful of large domains remain; the only 
dynamics beyond this is an extremely slow straightening 
of the interfaces. 

Surface roughness and domain size cannot provide such 
detailed information on the spatio-temporal dynamics. 
For a given structure, Wz,(f) saturates at the end of Stage 
1, and there is no discernible difference beyond this. The 
domain size grows as R(t) ~ fs at large times but also 
saturates between Stage 2 and 3, and gives no further 
information. 

We have also investigated the effect of adding zero- 
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FIG. 3: Measure for the KPZ model for different 5. L = 128. 



mean external noise to the CHE. We find that in the 
presence of noise, the dynamics in Stage 3 is faster than 
in the noise-free case, i.e., |<73(<5, t)\ is larger for the noisy 
dynamics. We also find that the saturation time of fi(S, t) 
is proportional to y^e. 

The KPZ equation is a paradigmatic model of non- 
equilibrium interfacial growth in the presence of lateral 
correlations 0, Q . It gives the local growth of the height 
profile h(x,t) at substrate position x and time t as 

^A = v v 2 h + \(vhf + V (x,t). (3) 

at 

Here v is the diffusion parameter which smoothens the 
interface and A the prefactor of the nonlinear term which 
tends to amplify large slopes. i)(x, t) is delta-correlated 
noise of zero mean which represents a random particle 
flux. RSOS is a discrete counterpart of the KPZ. Here, 
particles are added to a randomly chosen site i if and only 
if the addition ensures that all nearest-neighbor height 
differences \Ah\ < T, where T is some pre-determined 
positive integer. 

Roughness and correlation function analyses of the dy- 
namics of both KPZ and RSOS models show that their 
growth and roughness exponents take similar values in all 
dimensions, i.e., they are in the same universality class for 
large values of A. In the initial stages, the correlations 
spread across the entire lattice and the correlation length 
£ ~ ty , where z ~ 1.6. It has been established numeri- 
cally j2j) that the strong non- linear coupling limit of KPZ 
corresponds to the RSOS model in terms of the growth 
and the saturation exponents of the surface roughness. 

The behavior of fi(8, t) in time for the KPZ model 
(L = 128) is shown in Fig. The dynamics is quali- 
tatively similar to the behavior of the surface roughness, 
and can be divided into two distinct stages. In Stage 
1, the lateral correlations spread across the entire lat- 
tice and we find that fj,(8,t) decreases as t~ <Jl ( S ' t%> where 
cr\{8, t) is independent of 8 and t within numerical errors. 
The plots labelled by the dotted, dashed and solid lines 
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FIG. 4: The slopes vs. (5 for various A values compared with 
the RSOS slopes. 



in Fig. ^correspond to o\ (8, t) for three values of A of the 
KPZ model (A = 1, A = 2, and A = 25 respectively) dur- 
ing the first stage. [i(8, t) for the RSOS model exhibits a 
similar behavior during Stage 1. In Fig. 01 the plot with 
circles shows a\(8,t) versus 8 for the RSOS model. This 
corroborates the assertion that for the large A limit, KPZ 
corresponds to the RSOS models in terms of the growth 
exponents, since it is seen clearly that a\(8,t) for KPZ 
with A = 25 and for RSOS are within the errorbars. 

In Stage 2 of Fig|3 fi(S, t) saturates to a L-dependent 
value (given by L^W) where j(S) is the saturation ex- 
ponent. However, in this case, we find that -f(8) for 
KPZ and RSOS models, even at strong coupling, are 
very different . For example, 7(5) for the RSOS model is 
found to be 7(0.1) = -0.095 ± 0.004 while that for the 
KPZ model (A = 25) is 7(0.1) = -0.075 ± 0.007 when 
8 = 0.1. For 8 = 1.0, the saturation exponent for RSOS 
is 7(1.0) = -0.083 ± 0.002 while that for the KPZ is 
7(1.0) = -0.068 ±0.007. 



In summary, we have introduced a set of character- 
istics [1,(6, i) that can be used to study spatio-temporal 
dynamics of systems represented by a scalar field U(x, t). 
At a given instant, [1(8, t) are a set of inverse lengthscales 
associated with a pattern, and are defined in terms of the 
determinant of the Hessian. Large values of 8 emphasize 
regions with larger curvature of contour lines. 

The availability of a family of indices allows us to pro- 
vide a more comprehensive description than is possible 
from individual measures. For example, it was shown for 
the CHE that during Stage 2, relaxation of the distinct 
lengthscales associated with the structure occurs at dif- 
ferent rates. This is very different from Stage 1 where 
all such scales grow as fi(8,i) ~ t om . We are also able 
to identify instances where small domains of the pattern 
disappear by searching for peaks in fi(8, t) for large values 
of 8. 

Analysis of the KPZ and the RSOS models shows two 
stages in pattern evolution. For KPZ, [1(8, t) relaxes at a 
rate that depends on A during stage 1. The rate of decay 
<Ji(8,t) for large values of A (typically A > 15) is seen 
to be the same as for the RSOS model, thus reinforcing 
previous claims that both models belong to the same uni- 
versality class. However, our analysis provides an addi- 
tional piece of information, namely that all lengthscales 
[i^ 1 (S,t) associated with the spatio-temporal dynamics 
of these interfaces decay at the same rate. [i(8,t) sat- 
urates in stage 2 for both models, and the saturation 
value depends on the system size L as L 1 ^. Interest- 
ingly, the function j(8) for the two models is different; 
after saturation, the contours of the two models exhibit 
non-universal characteristics. It is thus possible to deter- 
mine which of these models better represents the growth 
of an experimental interface. 

The authors thank R. Rajesh and K. E. Bassler for a 
critical reading of the manuscript and for useful discus- 
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